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ABSTRACT 


The  introduction  of  gas  bubbles  into  a  liquid  decreases  the  average 
density,  and  thus  decreases  the  buoyant  force  on  a  floating  body.  This  thesis 
investigates  the  critical  average  density  required  to  sink  a  buoyant  body  in  water 
with  rising  bubbles.  A  volume  of  bubbly  water  is  created  in  a  clear  acrylic  tube  of 
inner  diameter  30  cm  and  height  60  cm,  that  is  closed  at  the  bottom  and  open  at 
the  top.  An  array  of  diffusers  at  the  bottom  produces  2  mm  diameter  bubbles 
distributed  uniformly  over  the  cross  section  of  the  tube.  A  1 0-cm  diameter  hollow 
steel  ball  whose  average  density  is  varied  from  0.70  to  0.99  g/cm3  is  employed 
as  the  buoyant  body.  A  theory  of  the  critical  density  for  sinking  is  developed,  and 
predicts  that  the  average  fluid  density  is  greater  than  the  ball  density  for  sinking. 
The  experimental  data,  which  include  a  quantitative  error  analysis,  agree  well 
with  the  theory  for  average  ball  densities  from  0.94  to  0.99  g/cm3,  but  show  a 
definite  trend  of  fluid  densities  that  are  smaller  than  those  predicted  for  0.70  to 
0.94  g/cm3. 
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I.  INTRODUCTION 


The  buoyant  force  on  a  body  that  is  partially  or  wholly  immersed  in  a  fluid 
in  a  gravitational  field  equals  the  weight  of  fluid  displaced  by  the  body 
(Archimedes’  Principle).  For  bodies  that  are  not  shaped  similar  to  a  canoe  or 
bowl,  which  can  displace  more  fluid  than  the  volume  of  the  material,  the  condition 
for  floating  is  simply  that  the  density  of  the  fluid  be  greater  than  the  average 
density  of  the  body.  Consider  such  a  body  floating  in  a  liquid.  The  introduction 
of  bubbles  results  in  a  fluid  whose  average  density  is  less  than  that  of  the  original 
liquid.  If  the  bubbles  are  uniformly  distributed  and  small  compared  to  the  size  of 
the  body,  we  might  expect  that  the  body  will  sink  when  the  average  density  of  the 
fluid  is  less  than  the  average  density  of  the  body.  However,  the  bubbles  may 
produce  upward  forces  on  the  body,  due  to  drag  produced  by  the  entrained  flows 
in  the  fluid  and  to  bubbles  sticking  to  the  body.  It  is  thus  not  clear  whether  the 
introduction  of  bubbles  can  make  a  floating  body  sink.  One  might  argue  that  the 
sinking  of  the  body  must  be  possible,  because  in  the  limiting  case  in  which  the 
gas  displaces  all  of  the  water,  any  object  more  dense  than  the  gas  would  sink. 

But  this  argument  is  faulty  because  it  neglects  a  possibly  very  large  upward  drag 
on  the  body,  which  could  prevent  the  body  from  falling. 

The  possible  sinking  due  to  bubbles  has  been  given  as  an  explanation  for 
the  sinking  of  some  ships,  especially  in  the  “Bermuda  Triangle”.  This  explanation 
is  based  on  the  existence  of  large  amounts  of  methane  gas  in  the  sediments  of 
the  region  (Rowe  and  Gettrust  1993).  This  gas  might  be  disturbed  and  released 
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in  sufficiently  large  quantities  such  that  the  density  of  the  water  is  temporarily 
reduced  to  levels  that  a  ship  cannot  displace  it  own  weight  of  water,  and  so  sinks 
(Mclver  1982).  Although  it  apparently  has  been  established  that  there  is  no 
greater  incident  of  sinkings  in  the  Bermuda  Triangle,  there  is  a  legitimate 
scientific  issue  regarding  the  sinking  of  a  floating  body  due  to  the  presence  of 
bubbles.  However,  we  have  not  found  any  literature  on  such  an  experiment. 
There  has  also  been  some  suggestion  that  the  deliberate  release  of  methane 
hydrate  deposits  by  underwater  vehicles  could  be  used  to  reduce  buoyancy  and 
thus  sink  targets.  This  is  referred  to  as  a  “Buoyancy  Bomb”  (Stumborg,  2000). 

Important  aspects  of  the  issue  are  the  uniformity  of  the  bubbles  and  the 
cross-sectional  area  over  which  the  bubbles  are  generated  compared  to  the 
cross-sectional  area  of  the  liquid.  Two  extreme  cases  are  that  (i)  the  bubbles  are 
generated  uniformly  over  the  cross-section  of  a  container,  and  (ii)  the  cross- 
sectional  area  of  the  liquid  is  much  greater  than  that  of  the  bubbles,  which  would 
occur  in  the  ocean,  for  example.  In  case  (i),  there  can  be  no  net  upward  flow  of 
liquid,  although  large-scale  nonuniformities  in  the  flow  could  conceivably  occur. 

In  the  absence  of  such  flows,  there  is  expected  to  be  little  if  any  drag  on  a  body 
due  to  the  liquid.  In  case  (ii),  large-scale  circulation  of  the  liquid  will  occur,  where 
the  liquid  rises  in  the  region  of  the  bubbles,  and  falls  outside  this  region.  This  is 
expected  to  produce  a  significant  upward  drag  on  a  body  in  the  region. 

The  purpose  of  this  thesis  is  to  describe  the  design  and  construction  of  an 
apparatus  that  will  allow  us  to  experimentally  investigate  case  (i),  and  to  compare 
the  data  to  theory.  We  have  observed  that  the  upward  drag  force  in  case  (ii)  is 
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substantial  and  we  recommend  that  these  results  be  investigated  in  future  work. 
We  show  that  a  body  floating  in  a  container  of  water  can  indeed  be  made  to  sink 
due  to  the  introduction  of  small  air  bubbles  uniformly  over  the  cross-sectional 
area  of  the  container.  This  offers  a  dramatic  lecture  demonstration,  which  is 
discussed  in  Appendix  A. 

The  organization  of  this  thesis  is  as  follows:  Chapter  II  addresses  the 
theory  determining  the  average  density  of  bubbly  water  at  which  a  body  is 
predicted  to  sink.  The  apparatus  that  we  have  constructed  for  this  purpose  is 
described  in  Chapter  III.  The  methods  of  calibration  for  our  system  are  provided 
in  Chapter  IV.  The  results  of  our  quantitative  investigations,  an  analysis  of  the 
system  uncertainties,  and  a  comparison  with  the  theory,  are  presented  in 
Chapter  V.  Concluding  remarks  and  suggestions  for  future  work  are  made  in 
Chapter  VI. 
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II.  THEORY 


We  require  a  method  of  determining  the  average  density  of  bubbly  water. 
Pure  water  of  volume  Vw  has  mass  pwVw,  where  pw  is  the  density  of  water.  If  air 
bubbles  of  total  volume  AV  are  introduced  into  the  water,  the  average  density  is 
then  p  =  pwVw/(Vw  +  AV),  where  we  neglect  the  mass  of  the  air  compared  to  the 
mass  of  the  water.  A  means  of  determining  this  density  is  to  measure  the  initial 
height  h0  of  the  water  level  in  the  container  with  no  bubbles,  and  the  change  in 
height  Ah  when  the  bubbles  are  turned  on  and  steady  state  has  been  reached.  If 
the  container  has  uniform  cross-sectional  area,  then  Ah/h0  =  AVA/W,  and  so  the 
average  density  of  the  bubbly  water  is: 


P 


Pw 

1  +  Ah/h0  ' 


(2.1) 


We  consider  the  critical  value  (2.1)  of  the  bubbly  water  density  below 
which  a  floating  ball  will  sink.  As  noted  in  the  Introduction  (Ch.  I),  if  bubbles  are 
produced  uniformly  over  the  cross-sectional  area  of  a  container,  there  will  be  no 
net  upward  flow  of  the  water  and  thus  no  upward  drag  on  the  ball.  The  critical 
value  should  then  ideally  equal  the  average  density  of  the  ball,  if  there  are  no 
large-scale  flows.  However,  above  the  ball  there  is  a  “shadow”  region  in  which 
there  are  no  bubbles,  whereas  this  critical  value  assumes  that  the  bubbles  are 
uniformly  distributed  throughout  the  entire  volume  of  liquid.  To  account  for  the 
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shadow,  we  assume  that  the  top  of  the  ball  is  tangent  to  the  surface  of  the  fluid 
when  the  ball  is  just  about  to  sink.  We  also  assume  that  the  shadow  has  uniform 
circular  cross  sectional  area  and  extends  vertically  from  the  equator  of  the  ball  to 
the  surface  of  the  fluid,  which  approximates  our  observations  of  the  system  (Ch. 
V.A).  The  volume  of  the  shadow  region  is  easily  shown  to  be  Vs  =  tiR3/3,  where 
R  is  the  radius  of  the  ball.  The  absence  of  bubbles  above  the  ball  causes  the 
pressure  on  the  top  of  the  ball  to  be  greater  than  it  would  be  if  bubbles  were 
present.  This  force  is  the  weight  of  the  additional  water  above  the  ball  due  to  the 
lack  of  bubbles,  so  there  is  a  downward  force  (pw  -  p)gVs  on  the  ball. 

Substitution  of  the  above  expression  for  Vs  yields  the  downward  force  on  the  ball 
due  to  the  shadow  region  as 

Fs  =  §(A*  -  p)9P?  •  (2.2) 

Because  the  ball  alters  the  flow  of  the  bubbles,  the  bubble  density  near 
the  bottom  of  the  ball  and  the  cylindrical  boundary  of  the  shadow  region  is 
conceivably  greater  than  that  farther  from  the  ball,  which  may  cause  a  significant 
additional  decrease  of  the  buoyancy.  However,  our  observations  (Ch.  V.A) 
reveal  that  the  bubble  density  is  approximately  uniform  throughout  the  fluid.  The 
bubbles  that  would  have  been  in  the  shadow  region  are  thus  effectively  uniformly 
distributed  throughout  the  rest  of  the  water,  which  reduces  the  buoyancy  by  a 
negligible  amount  in  the  typical  case  where  the  volume  of  the  water  is  much 
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greater  than  the  volume  of  the  ball. 

The  condition  for  neutral  buoyancy  of  the  ball  is  W  +  Fs  =  FB,  where  the 
downward  gravitational  force  is  W  =  mg  and  the  upward  buoyant  force  for  a 
uniform  fluid  is  FB  =  pVg,  where  m  is  the  mass  of  the  ball  and  V  is  its  volume. 
Substituting  the  expression  (2.2)  of  the  shadow  force  Fs,  and  solving  for  the 
average  fluid  density  p,  yields  the  value  below  which  the  ball  is  predicted  to  sink: 


m  +  pwnR3/3 
P  V  +  7tR3/ 3 


(2.3) 


This  expression  correctly  reduces  to  pw  =  m/V  in  the  case  of  no  bubbles  (p 
=  pw).  We  do  not  set  V  =  4tcR3/3  in  Equation  (2.3)  because  the  ball  used  in  the 
experiment  (Ch.  V.A)  has  a  small  protrusion  as  part  of  the  filling  hole.  The 
protrusion  points  downward  in  the  experiment,  and  causes  the  volume  V  of  the 
ball  to  be  0.6%  greater  than  4tcR3/3.  The  ball  is  also  slightly  oblate  by  0.25%. 

We  minimize  the  error  due  to  the  oblateness  by  substituting  the  average  radius 
for  R  in  Equation  (2.3). 

It  should  be  noted  that  the  above  theory  predicts  that  the  equilibrium  state 
(in  which  the  top  of  the  ball  is  tangent  to  the  surface  of  the  fluid)  is  unstable.  This 
occurs  because  a  small  downward  displacement  of  the  ball  causes  the 
downward  force  due  to  the  shadow  to  increase,  but  does  not  alter  the  buoyant 
force  for  a  uniform  fluid. 
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Figure  3.1 .  Main  components  of  the  apparatus:  a  cylindrical  acrylic  chamber,  and  an  array  of  1 5 
bubble  diffusers.  Also  shown  is  the  buoyant  body,  which  was  employed  throughout  the  thesis. 
The  body  has  been  sunk  due  to  bubbles.  To  avoid  distortion,  this  picture  was  taken  with  the 
chamber  submerged  in  a  large  reservoir  of  water. 
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A.  THE  BUBBLY  WATER  CHAMBER 


The  bubbly  water  environment  is  created  within  an  acrylic  tube  that  is 
closed  at  the  bottom  and  open  at  the  top.  The  dimensions  of  the  tube  are  30.6 
cm  for  inner  diameter  and  60  cm  in  height  with  a  wall  thickness  of  3.175  mm. 

The  water  depth  is  maintained  at  approximately  40  cm  to  allow  for  a  nearly  20  cm 
increase  in  height  resulting  from  the  introduction  of  the  bubbles. 

The  chamber’s  cylindrical  shape  was  selected  to  minimize  the  wall  surface 
area  for  a  given  volume  and  thus  minimize  the  surface  drag  of  the  bubbles  on  the 
walls  of  the  chamber.  The  cylindrical  shape  also  provided  no  corners  that  might 
needlessly  complicate  the  bubble  flow  pattern.  The  diameter  of  the  chamber  was 
chosen  to  be  as  large  as  conveniently  possible  compared  to  the  diameter  of  the 
buoyant  body.  The  chamber  height  allows  for  a  uniform  mixing  of  the  bubbles 
ensuring  an  equal  distribution  of  the  small  air  bubbles  over  the  cross-sectional 
area  of  the  chamber.  The  clear  acrylic  material  facilitates  viewing  and 
photography  of  the  bubbles  and  the  buoyant  body. 

Holes  of  diameter  0.7  cm  were  drilled  3  cm  from  the  bottom  of  the 
chamber  to  allow  15  bubble-generating  diffusers  to  penetrate  the  cylinder.  The 
holes  were  placed  such  that  the  diffusers  were  evenly  spaced  and  parallel  to 
each  other  (Fig.  3.2).  To  determine  where  to  place  the  holes,  we  required  the  arc 
length  to  the  center  of  diffuser  from  an  initial  point  on  the  tube.  This  was  found  to 
be: 
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ArcLength(N)  =  RTCos  1 

where  tfr  =  tube  radius  (OD)  =  153  mm,  tT  =tube  thickness  =  3.175  mm,  dT  = 
distance  of  the  1st  diffuser  from  the  chamber  wall  =  4  mm,  Dd  =  diffuser  diameter 
=  5.6  mm,  N  =  diffuser  number  =  1-15,  and  Sd  =  distance  from  center  to  center 
of  diffusers  =  20  mm. 


Figure  3.2.  Top  view  of  the  bubbly  water  chamber  and  diffuser  arrangement. 

A  6  mm  thick  acrylic  plate,  33  cm  square,  was  glued  to  the  bottom  of  the 
tube.  Once  the  diffusers  were  inserted  into  the  holes,  the  plate  effectively 
created  a  chamber  that  was  both  water  and  air  tight  below  the  water  line.  This 
proved  to  be  an  important  feature  during  the  acquisition  of  the  experimental  data 
(Ch.  V.A). 


WT 


(3.1) 
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B.  VOLUMETRIC  BUBBLE  GENERATOR 


Several  bubble  making  devices,  including  perforated  coiled  flexible  tubing, 
“frits”  and  perforated  acrylic  plates,  were  considered  prior  to  the  final  selection 
being  made.  The  three  primary  criteria  for  the  diffuser  selection  were:  the 
bubbles  must  be  much  smaller  than  the  buoyant  body  yet  large  enough  so  that 
direct  absorption  of  the  nitrogen  gas  into  the  fluid  is  inhibited;  there  must  be  a 
constant  and  diffused  flow  of  bubbles  to  ensure  a  uniform  distribution  over  the 
cross-sectional  area  of  the  chamber,  reducing  possible  fluid  turbulence;  and  the 
diffusers  must  be  sufficiently  robust  to  physically  withstand  a  large  volumetric 
flow  of  nitrogen  over  extended  periods.  The  device  that  provided  the  most 
uniform  bubble  field,  with  small  diameter  bubbles  and  little  turbulence,  was  a  set 
of  aquarium  diffusers  (Regent  PL-T714).  The  diffusers  are  hollow  porous  tubes 
capped  at  one  and  with  an  air  inlet  at  the  other.  The  dimensions  of  the  diffusers 
are  36  cm  in  length,  5.6  mm  in  diameter  and  a  wall  thickness  of  .75  mm. 

A  visual  examination  of  individual  Regent  diffusers  showed  a  wide 
variation  of  performance  between  diffusers,  both  in  pattern  flow  and  in  bubble 
size.  A  Canon  Optura  digital  video  camera  was  employed  to  examine  the  flow  of 
25  individual  diffusers,  for  longitudinal  uniformity  and  bubble  size  uniformity. 
Fifteen  diffusers  were  selected  from  this  lot  and  strategically  positioned  to  give 
the  most  uniform  cross-sectional  flow  pattern  possible.  The  diffusers  were 
inserted  through  the  chamber  holes  (Sec.  A).  The  portion  of  the  diffusers  that 
extended  outside  of  the  chamber  was  made  water  and  air  tight  through  the  use  of 
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heat  shrink  and  silicon  sealant. 


C.  NITROGEN  GAS  FEED  SYSTEM 

A  schematic  of  the  feed  system  is  shown  in  Figure  3.3.  Accurate  control 
of  the  volume  flow  rate  to  the  diffusers  was  important.  Overall  gas  flow  rate  was 
controlled  with  a  single  stage  pressure  regulator  (Matheson  Gas  Products  Model 
1L-580)  in  conjunction  with  an  inline  flow  meter.  Tygon  tubing  (0.3125  in  ID) 
connected  the  nitrogen  source  to  individual  diffuser  control  valves  (Copper  - 
9075656)  through  a  manifold  (locally  fabricated).  Balancing  of  the  diffuser 
control  valves  was  found  to  be  critical  in  the  stability  of  the  generation  of  bubbles. 
Any  small  deviation  from  a  perfect  balance  set  up  substantial  turbulent  flow  in  the 
chamber  (Ch.  IV.A.5).  It  is  recommended  that  multi-turn  needle  valves  be 
utilized  in  the  future  to  replace  the  inexpensive  gate  valves  that  we  employed. 


Pressure  Gauges 


Nitrogen  Cylinder  Control  Valves  Diffusers 

Figure  3.3.  Schematic  of  the  nitrogen  gas  feed  system.  There  are  15  diffusers. 
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IV.  METHODS  AND  SYSTEM  CALIBRATION 


A.  METHODS 

As  explained  in  Chapter  II,  if  air  bubbles  of  total  volume  AV  are  introduced 
into  the  water,  the  average  density  is  p  =  pwVw/(Vw  +  AV).  The  ability  to 
accurately  determine  the  density  of  the  bubbly  water  was  critical  to  our 
experiment.  A  method  that  consistently  provided  a  means  to  accurately 
determine  the  density  proved  to  be  a  challenge  and  required  several  evolutions. 
Five  different  methods  to  determine  the  average  density  of  the  bubbly  water  were 
attempted,  one  with  very  good  outcomes,  others  with  limited  or  negligible  results. 

1 .  Bubble  Velocity  Method 

The  volume  of  gas  that  passes  through  the  flowmeter  in  a  certain  amount 
of  time  T  is  the  product  of  the  volume  flow  rate  Q  and  the  time.  When 
considering  the  volume  of  gas  AV  contained  in  the  bubbles  in  the  water,  the 
relevant  time  is  that  required  for  the  bubbles  to  rise  from  the  bubble  generators  to 
the  surface.  To  determine  AV,  we  note  that,  as  a  result  of  viscous  losses  in  the 
tubes  and  manifold,  the  pressure  pt  at  the  flowmeter  is  greater  than  the  pressure 
at  the  bottom,  which  is  approximately  atmospheric  pressure  p0.  We  make  the 
assumption  that  the  air  remains  at  constant  temperature  as  it  flows,  so  that  the 
product  of  pressure  and  volume  is  constant  (Boyle’s  law).  This  yields  ptQT  = 
p0AV.  Solving  for  AV  yields: 
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pw 


(4.1) 


1  +  p,QT/p0Vw  ' 

Equation  (4.1)  allows  the  average  density  of  bubbly  water  to  be  calculated  from 
parameters  that  can  be  measured. 

A  digital  video  camera  (Canon  Optura)  was  employed  to  capture  the  time- 
distance  sequence  of  a  typical  bubble  of  2  mm  diameter.  A  meter  stick 
photographed  in  the  plane  of  the  bubble  provided  the  distance  scale  factor  for  the 
evaluation  of  displacement  and  size.  The  average  vertical  speed  was  found  to 
vary  from  25.2  to  34.3  cm/s  depending  on  the  flow  rate.  This  is  consistent  with 
Morton  (1953)  who  concluded  that  a  single  2  mm  bubble  reached  a  velocity  of  24 
cm/s  in  water,  and  that  the  presence  of  wake  in  the  liquid  (caused  by  other 
bubbles  in  close  proximity)  results  in  higher  velocity  of  rise  of  the  bubble.  Our 
bubble  velocity  method  proved  to  be  accurate  only  to  a  volume  flow  rate  of  3 
l/min,  at  which  the  larger  quantity  of  bubbles  made  it  very  difficult  to  discern  one 
bubble  from  another.  The  3  l/min  flow  rate  represented  only  about  2%  of  the 
expected  range  of  our  experiment.  The  problem  is  that  we  cannot  assume  that 
the  bubble  velocity  is  independent  of  flow  rate,  because  the  bubble  size  may 
change,  and  the  bubbles  may  interact.  Therefore  other  methods  needed  to  be 
explored. 

2.  Trailing  Edge  Method 

An  alternative  approach  to  the  above  method  was  attempted.  We 
measured  the  rise  time  T  by  stopping  the  airflow  and  timing  the  trailing  edge  of 
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the  bubbles  as  it  passed  from  a  lower  plane  to  the  surface.  Although  the  bubble 
diameters  averaged  2  mm,  they  in  fact  ranged  from  0.5  mm  to  3  mm  and 
consequently  the  different  velocities  made  a  discrete  trailing  edge  too  difficult  to 
distinguish. 

3.  Acoustical  Method 

The  speed  of  sound  in  a  fluid  is  given  by: 


c 


2 


B^=  1 

P  PKad  ’ 


(4.2) 


where  Bad  is  the  adiabatic  bulk  modulus,  p  is  the  mass  density  and  Kad  =  B~\d  is 

the  adiabatic  compressibility  of  water  (Crocker,  1998).  If  an  amount  of 
undissolved  nitrogen  gas  is  uniformly  mixed  with  water,  the  mixture  may  be 
treated  as  a  homogenous  fluid  for  sufficiently  large  wavelength.  Equation  (4.2) 
may  then  be  applied,  where  the  density  and  compressibility  are  the  volume- 
weighted  averages  of  the  gas  and  the  liquid  densities  and  compressibilities: 


p  =  +  (l  -  X)P,  , 

(4.3) 

Kad  =  XKg+{\-X)K,. 

(4.4) 

Here  the  subscripts  g  and  /  refer  to  gas  and  liquid,  respectively,  and*  is  the 
gas  volume  fraction.  Substituting  Equations  (4.3)  and  (4.4)  into  Equation  (4.2) 
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and  solving  for  x  yields: 


X  =  - 


c2{Kl^-2K9 


KiPgPl  +  Kfp2g)  +  A(Kg  -Kj){Pg -Pl)-C{KgPi  +  Kj{Pg  -2p/)) 

2c{Kg-Kl){Pg-Pj) 


(4.5) 


Equation  (4.5)  allows  x  to  be  determined  by  measuring  c. 

A  point  source  transducer  (USN  DT-574/BQQ-6)  was  driven  by  a  3-pulse 
burst  from  a  1 5  MHz  Function/Arbitrary  Waveform  generator  (HP  331 20A) 
through  a  power  amplifier  (HP  467A)  at  the  transducer’s  resonant  frequency  of 
26.65  kHz.  Two  pickup  hydrophones  (EDO  Model  6600)  were  positioned  in  line 
with  the  source  beam.  The  receive  signals  were  recorded  on  an  oscilloscope 
(Agilent  Infiniium  5481 0A)  through  low  noise  pre-amplifiers/filters  (Stanford 
Research  Systems  Model  SR  560).  Knowing  the  distance  between  the 
hydrophones  and  measuring  the  time  between  the  two  received  waveforms  allow 
us  to  easily  determine  the  speed  of  sound  in  water  without  bubbles.  However, 
when  bubbles  were  introduced  into  the  water,  a  frequency  spectrum  shows  a 
loss  greater  than  70  dB  in  the  received  signals  and  as  such  the  transmitted  signal 
became  ‘buried’  in  the  ambient  noise.  Increasing  the  source  level  and  further 
amplifying  the  received  signal  provided  no  further  assistance  in  raising  the 
receive  signal  out  of  the  noise. 

We  ruled  out  attenuation  of  the  signal  due  to  resonance  frequency  which 
is: 
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(4.6) 


= j_  fiK 

2 na\  p 


with  a  =  bubble  radius,  p0  =  density  for  water,  y  =  specific  heat  constant  and 
P6  =  P0  (1  +  z/10)  for  z  in  meters  (Kinsler  et  al.  2000).  For  our  average  2  mm 

diameter  bubble,  at  a  depth  of  10  cm,  the  resonant  frequency  is  3.28  kHz,  which 
is  much  less  than  the  27  kHz  frequency  of  the  sound.  The  lack  of  signal  at  the 
pickup  hydrophones  could  be  due  to  ineffective  coupling  of  the  driving  transducer 
into  the  medium,  or  to  the  scattering  of  the  sound  waves  off  the  bubbles  or 
turbulence  caused  by  the  bubbles.  This  technique  was  thus  ruled  out  as  an 
inaccurate  means  of  measuring  the  density  of  our  bubbly  water. 

4.  Dielectric  Material  Method 

By  measuring  the  capacitance  C  of  a  parallel  plate  capacitor  with  the 
space  between  the  plates  having  area  A  filled  with  bubbly  water  and  the 
capacitance  Cw  with  water  between  the  plates,  we  can  find  the  average  density 
of  the  bubbly  water  in  the  chamber.  From  the  relationship  of  capacitance 

K  £  A 

C  =  — — ,  where  e  is  the  permittivity  of  free  space,  Equation  (4.3),  and 
d 

Kd  -  XKg  +  (1  -  x)k ) ,  we  find: 

C 

P  =  Pw-^r,  (4.7) 
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where  Kg  and  k,  is  the  dielectric  constant  of  the  gas  and  liquid  respectively,  and 
with  the  understanding  that  k ,  »  k9  ,  and  p,  »  pg . 

To  acquire  an  accurate  average  density  of  the  bubbly  water,  it  was 
estimated  that  the  plates  would  have  to  be  at  least  15  cm  apart.  With  this  plate 
separation,  the  largest  plate  area  that  would  fit  into  the  chamber  was  calculated 
to  be  600  cm2.  The  theoretical  capacitance  of  water  for  a  parallel  plate  capacitor 
with  the  above  parameters  is  273  pF. 

A  preliminary  test  to  determine  the  practicality  of  this  method  was 
conducted.  Two  unetched,  copper  clad  printed  circuit  boards,  with  dimensions 
20  cm  x  30  cm,  were  placed  in  parallel  and  25  cm  apart  on  the  outside  walls  of  a 
glass  aquarium.  The  aquarium  was  filled  with  water.  The  copper  plates  were 
connected  to  an  impedance  analyzer  (HP  41 94A)  using  a  two-terminal 
measurement  configuration.  The  results  of  our  test  were  inconsistent  with  theory 
and  did  not  vary  greatly  when  bubbles  were  introduced  to  the  water.  We 
attributed  this  to  two  assumptions  in  Equation  (4.7):  the  space  between  the 
plates  is  filled  with  a  homogenous  dielectric;  and  the  plate  separation  is  small 
compared  with  the  dimensions  of  the  plate.  The  latter  would  allow  the  fringing 
field  of  the  plates  to  be  ignored.  Both  of  these  assumptions  were  violated  in  our 
test  set  up  and  it  was  predicted  that  similar  results  could  be  expected  if  employed 
in  the  chamber.  This  technique  was  also  ruled  out  as  an  accurate  means  of 
measuring  the  density  of  our  bubbly  water. 
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5.  Change  in  Height  Method 

A  standard  means  of  determining  the  density  of  bubbly  water  is  to 
measure  the  initial  height  h  of  the  water  level  in  the  chamber  with  no  bubbles, 
and  the  change  in  height  Ah  when  the  bubbles  are  turned  on  and  steady  state 
has  been  reached.  If  the  container  has  uniform  cross-sectional  area  the  average 
density  of  the  bubbly  water  can  be  determined  from  Equation  (2.1 ). 

We  found  that  we  were  able  to  accurately  measure  the  change  in  height 
only  after  the  diffusers  were  in  perfect  balance;  that  is,  the  flow  rate  per  unit 
length  was  identical  over  each  diffuser.  When  slightly  out  of  balance  the  bubble 
flow  became  very  turbulent,  eddies  were  observed  within  the  volume  of  bubbly 
water  and  the  surface  took  on  a  boiling  like  appearance.  We  were  then  unable  to 
measure  the  change  in  height  with  any  sort  of  precision. 

The  nitrogen  source  originally  fed  all  the  diffusers  from  the  same  direction. 
Visually  it  appeared  that  there  was  a  greater  flow  of  bubbles  stemming  from  the 
source  end  than  from  the  opposite  end.  By  alternating  the  source  direction  for 
every  other  diffuser  the  turbulence  was  reduced  significantly.  Further  reduction 
was  accomplished  by  meticulously  adjusting  each  of  the  15  control  valves  such 
that  the  bubble  flow  from  each  diffuser  was  visually  balanced.  The  bubble  flow 
became  uniform  throughout,  the  eddies  completely  disappeared,  and  the  surface 
became  level. 
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B.  SYSTEM  CALIBRATION 


The  ability  to  accurately  measure  the  change  in  height  (Sec.  A.5)  allowed 
us  to  use  Equation  (2.1)  to  determine  the  density  of  the  bubbly  water.  To  simplify 
this  procedure  we  calibrated  the  change  in  height  Ah  as  a  function  of  the  total 
gas  flow  rate  Q.  Flow  rates  ranging  from  0  to  134  l/min,  at  15  different  values, 
were  set  at  the  nitrogen  source,  and  the  corresponding  changes  in  height  were 
recorded.  Three  different  flowmeters  were  required  to  cover  the  range  of 
interest.  Two  150-mm  direct-reading  flowmeters  (Gilmont  Instruments)  covered 
the  lower  range  0-16  l/min  (GE717)  and  16  -  44  l/min  (GE817),  and  a  digital 
flow  sensor  (McMillian  1 00-1 8  Flo-Sen)  was  employed  at  the  higher  flow  rates. 
The  results  of  the  system  calibration  are  shown  in  Figure  4.1 .  A  3rd  order 
polynomial  was  used  in  a  least  squares  fit  through  the  points,  resulting  in: 

A/7  =  a  +  b1Q+b2Q2  +  Jb3Q3,  (4.8) 


where  a  =  -0.171507,/}  =  0.419662, b2  =3.330407x1 0-4,  and 

b3  =  3.297877 x  1 0"5 .  The  flow  rate  Q  is  measured  in  l/min  and  change  in  height 
Ah  is  measured  in  mm. 
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V.  EXPERIMENT 


Our  objective  in  the  experiment  is  to  make  a  buoyant  body  sink  by  adding 
bubbles,  and  to  gather  the  data  needed  to  use  Equation  (2.1)  to  calculate  the 
density  of  the  fluid  when  this  occurs.  We  conduct  an  error  analysis  on  the 
system  and  finally  discuss  the  results  of  our  findings. 

A.  EXPERIMENT 

We  required  a  buoyant  body  that  was  much  larger  relative  to  the  bubbles 
so  that  we  could  treat  the  effects  of  the  bubbly  water  on  the  object  as  those  of  a 
homogenous  fluid.  We  selected  a  commercially  available  “density  ball”,  which  is 
a  hollow,  stainless  steel,  101.3  mm  diameter  sphere  with  a  removable  plug 
(Cenco  Scientific  Co.  cat.  #76595N).  Liquid  can  be  added  or  removed  to  vary 
the  ball’s  mass  and  thus  its  average  density.  We  employed  mineral  oil  in  order  to 
prevent  rusting. 

We  obtain  an  accurate  value  of  the  volume  of  the  ball  by  a  variation  of 
Archimedes’  method  (weighing  a  body  in  and  out  of  water).  Sufficient  liquid  is 
added  to  the  ball  such  that  it  will  sink  in  water.  We  then  place  a  container  of 
water  on  a  digital  balance  (AND  Electronic  Balance  FX-2000),  and  “zero”  the 
balance.  The  ball  is  suspended  by  string  from  a  support  and  is  totally  immersed 
in  the  water  without  touching  the  container.  We  record  the  reading  m  of  the 
scale.  The  buoyant  force  mg  equals  the  weight  pwVbaii9  of  the  displaced  water, 
so  the  volume  of  the  ball  is  Vbaii  =  m/pw  =  556.73  cm3.  This  method  has  two 
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advantages  over  Archimedes’  method:  arranging  suspension  of  a  body  from  an 
elevated  balance  is  avoided,  and  the  buoyancy  due  to  air  is  irrelevant.  The  mass 
of  the  empty  ball  is  229.61  g,  so  the  minimum  average  density  of  the  ball  is 
0.4124  g/cm3 

For  the  main  part  of  the  experiment,  we  determine  the  average  density  of 
the  bubbly  water  required  to  sink  the  ball.  We  select  an  amount  of  liquid  in  the 
ball  such  that  it  floats  in  the  water.  We  place  the  ball  in  the  chamber,  slowly 
increase  the  flow  rate  Q,  and  record  the  minimum  value  for  which  the  ball  sinks 
below  the  surface.  We  then  compare  this  to  the  ball’s  average  density,  which  is 
determined  as  the  total  mass  (including  liquid  inside)  divided  by  the  volume  of  the 
ball.  The  entire  process  is  then  repeated  for  a  different  value  of  the  average 
density  of  the  ball. 

In  our  initial  experiments,  we  observed  two  effects  that  cause  deviations 
from  the  measured  values  of  the  requisite  airflow  rates,  but  which  could  be 
eliminated.  First,  bubbles  sticking  to  the  lower  half  of  the  ball  effectively  lower 
the  ball’s  average  density,  which  make  the  ball  neutrally  buoyant  at  a  higher 
airflow  rate  than  expected.  Second,  when  the  ball  is  floating  at  the  surface,  there 
is  a  meniscus  that  forms  around  the  ball  at  the  ball-water  interface,  so  surface 
tension  is  acting  upward  on  the  ball.  This  force  must  be  overcome  with  a  higher 
airflow  rate  when  sinking  the  ball.  The  number  of  bubbles  sticking  to  the  surface 
of  the  ball  was  reduced  to  nearly  zero  by  machine  polishing  the  ball  and 
thoroughly  cleaning  it  with  isopropyl  alcohol  prior  to  each  trial.  Considerable  care 
was  also  taken  when  handling  the  ball  to  prevent  any  fingerprints  on  the  surface. 
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We  substantially  reduced  the  surface  tension  by  applying  to  the  outer  ball  surface 
a  very  thin  coat  of  Photo-Flo  200,  which  is  a  surfactant  used  in  film  developing. 
Another  effect,  the  absence  of  bubbles  above  the  ball,  could  not  be  eliminated, 
so  we  incorporated  estimations  of  this  effect  into  the  theory  (Ch.  II). 

The  criterion  we  set  for  determining  the  sinking  of  the  ball  was  that  it 
remain  below  the  surface  of  the  water  for  roughly  70%  of  the  time  over  a  30 
second  period.  This  benchmark  was  necessary  to  define  the  sinking  of  the  ball 
because  as  the  fluid  density  approached  the  density  of  the  ball,  the  ball  tended  to 
“loiter”  in  a  volume  of  bubbly  water  just  below  the  surface  of  the  water  rather  than 
sinking  directly  to  the  bottom.  Occasionally,  during  this  loitering  phase,  the  ball 
would  subtly  break  the  surface  of  the  water. 

Figure  (5.1)  shows  the  results  of  an  experiment  to  determine  the  average 
fluid  density  at  which  the  ball  sinks.  The  fluid  density  at  which  the  ball  just  sinks 
below  the  surface  is  significantly  greater  than  the  density  of  the  ball,  which  is 
represented  by  the  identity  line  in  the  graph.  A  3rd  order  polynomial  was  used  in 
a  least  squares  fit  through  the  points,  resulting  in: 

—  =  c  +  d,Q+  d2C?  +  d3Cf ,  (5.1) 

Pw 

where  c  =  0.999826,  d,  =  -8.694053x1  O'4 ,  d2  =- 4.97291 1x1  O’6 , 
d3  =  -2.06743 xIO"8,  and  where  Q  is  measured  in  l/min. 
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Figure  5.1 .  The  fluid  density  at  which  the  ball  just  sinks  below  the  surface  is  significantly  more 
than  the  density  of  the  ball,  which  is  represented  by  the  identity  line  in  the  graph. 

Several  special  observations  were  noted  during  our  experiment.  First,  it 
appeared  to  all  observers  that  the  volume  fraction  of  bubbles  was  much  greater 
than  what  it  actually  was,  judging  solely  by  the  vast  number  of  bubbles  in  the 
chamber.  Only  after  determining  the  actual  density  using  the  change  in  height 
method  (Ch.  IV.A.5)  was  an  observer  convinced  that  his  or  her  visual  intuition 
was  incorrect.  Second,  as  the  flow  rate  increased  above  5  l/min  the  bubble 
velocity  decreased.  This  is  contrary  to  Morton’s  (1953)  observation  that  bubble 
velocity  increases  in  the  presence  of  wake.  Solving  for  T  in  Equation  (4.1),  the 
bubble  velocity  v  =  (h0  +  Ah)/T  results  in  our  finding  shown  in  Figure  (5.2).  Third, 
when  the  top  of  the  ball  is  tangent  to  the  surface  of  the  fluid,  there  exists  a 
shadow  above  the  ball  in  which  there  are  no  bubbles.  The  shadow  has  a  uniform 


28 


circular  cross-sectional  area  that  extends  vertically  from  the  equator  of  the  ball  to 
the  surface  of  the  fluid.  Once  the  ball  is  sunk  the  “plume”  bulges  outward. 


Figure  5.2.  Above  5  l/min  the  bubble  velocity  decreased  which  is  contrary  to  Morton’s 
observation  that  bubble  velocity  increases  in  the  presence  of  wake. 


B.  ERROR  ANALYSIS 

There  are  a  number  of  sources  of  error  in  the  determination  of  the  density 
at  which  the  ball  sinks.  Uncertainties  exist  in  the  flowmeters  outputs,  the 
measurement  in  the  change  of  the  surface  height,  and  the  flow  rate  at  which  the 
ball  sinks  below  the  surface.  In  this  analysis,  all  uncertainties  are  worst-case 
tolerances. 

We  first  consider  the  system  calibration  (Ch.  IV.B),  for  which  there  are 
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three  independent  contributions  to  the  uncertainty:  the  repeatability  of  the 
flowmeter,  the  resolution  of  the  flowmeter,  and  the  uncertainty  8h  in  measuring 
the  surface  height  after  the  introduction  of  bubbles.  Both  the  repeatability  and 
resolution  are  taken  directly  from  the  flowmeter  manufacturers’  specifications 
(App.  B).  The  Gilmont  flowmeters  are  flow  rate  dependent.  The  McMillian 
flowmeter  is  not.  Although  we  found  the  surface  to  be  level  during  system 
calibration,  there  is  some  uncertainty  when  measuring  the  surface  height  h.  This 
we  determined  to  be  8h  =  ±0.5  mm  regardless  of  flowrate  Q.  The  result  of  our 
system  calibration  error  analysis  is  shown  as  error  bars  in  Figure  5.3. 


Figure  5.3.  Calibration  curve  for  the  change  in  height  Ah  as  a  function  of  the  flow  rate  Q. 
Independent  uncertainties  exist  in  Ah  directly  and  in  Q.  The  latter  indirectly  causes  an  uncertainty 
in  Ah,  which  is  negligible  on  the  scale  of  the  figure. 


We  next  consider  the  experiment  (Sec.  A),  in  which  we  determine  the 
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specific  gravity  of  the  fluid  at  which  the  ball  sinks.  There  is  one  independent 
contribution  to  the  uncertainty:  the  repeatability  of  the  experiment.  Seven 
separate  trials  were  conducted  for  a  given  specific  gravity  of  the  ball.  The  flow 
rate  for  which  the  ball  sinks  below  the  surface  was  recorded  each  time.  The 
uncertainty  in  the  experiment  repeatability  was  calculated  as  half  of  the 
difference  between  the  maximum  and  minimum  flow  rates  recorded. 

The  overall  uncertainty  in  the  flow  rate  8Q  is  the  quadrature  sum  of  the 
flowmeter  repeatability,  the  flowmeter  resolution,  and  the  experiment 
repeatability.  The  error  bars  in  Figure  5.4  show  the  uncertainty  of  the  specific 
gravity  of  the  fluid,  which  is  a  direct  function  of  the  uncertainty  in  the  change  in 
height  5(Ah),  and  is  determined  by: 


8  (Ah)  =  yj(8h)2  +  (df(A  h)ldOf(8Q)z,  (5.2) 

where  f(Ah)  is  a  3rd  order  polynomial  least  squares  fit  (Eq.  4.8),  and  where  the 
uncertainty  of  the  specific  gravity  is: 

Jo  =  1  S(Ah) .  (5.3) 

Pw  (1  +  Ah/hg)  h0 

Definitions,  specifications  of  flowmeters,  and  detailed  results  of 
uncertainties  are  provided  in  Appendix  B. 
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Figure  5.4.  Uncertainty  of  the  specific  gravity  of  the  fluid  is  a  direct  function  of  the  uncertainty  in 
the  change  in  height  8(Ah). 


C.  SUMMARY  OF  FINDINGS 

The  results  of  our  experiment  clearly  shows  that  the  fluid  density  at  which 
the  ball  just  sinks  below  the  surface  is  significantly  greater  than  the  density  of  the 
ball,  which  is  represented  by  the  identity  line  in  the  graph  (Fig.  5.4).  The 
experimental  data,  which  include  the  quantitative  error  analysis,  agree  well  with 
the  theory  for  average  ball  densities  from  0.94  to  0.99  g/cm3,  but  show  a  definite 
trend  of  fluid  densities  that  are  smaller  than  those  predicted  for  0.70  to  0.93 
g/cm3  (Fig.  5.5  and  5.6). 
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Specific  Gravity  of  Ball 


Figure  5.5.  The  fluid  density  at  which  the  ball  just  sinks  below  the  surface  is  significantly  greater 
than  the  density  of  the  ball,  which  is  represented  by  the  identity  line  in  the  graph.  The 
experimental  data  agree  well  with  the  theory  for  specific  gravity  of  the  ball  from  0.94  to  0.99,  but 
show  a  definite  trend  of  fluid  densities  that  are  smaller  than  those  predicted  for  0.70  to  0.93. 


Figure  5.6.  An  expanded  view  of  Fig  (5.5),  showing  the  region  of  transition  from  dilute  Newtonian 
fluid  behavior  to  possibly  non-Newtonian  fluid. 
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We  consider  two  possibilities  for  this  abrupt  transition  between  0.94  -  0.93 
g/cm3.  First  we  suspect  that  there  may  exist  large-scale  nonuniformities  in  the 
flow  as  discussed  in  Chapter  I.  To  investigate,  we  placed  in  the  chamber  a 
spherical  body  that  is  much  smaller  than  the  chamber  diameter  and  with  a 
density  of  0.91  g/cm3.  We  slowly  increased  the  flow  rate  Q  and  stopped  a 
number  of  times  prior  to  the  body  actually  sinking.  During  these  stopped  times 
we  carefully  observed  the  response  of  the  body  as  it  explored  the  cross-sectional 
area  of  the  chamber.  No  changes  in  the  water  depth  of  the  body  were  witnessed 
thus  we  conclude  that  there  are  no  large-scale  nonuniformities  in  the  flow. 

Second  we  speculate  that  the  abrupt  transition  may  be  linked  to  a  change  of 
flowmeters  near  the  flow  rate  region  of  interest.  However,  a  third  flowmeter 
(McMillian  100-13  Flow-Sen)  covering  a  range  of  20-100  l/min  was  employed  to 
verify  previous  data.  The  results  were  nearly  identical  thus  we  deduce  that  the 
change  of  flowmeters  was  not  the  cause. 

The  abrupt  transition  may  be  attributed  to  a  change  from  Newtonian  to 
non-Newtonian  behavior  of  the  fluid  (Wallis,  1969).  At  low  flow  rates,  we  can 
treat  the  fluid  as  simply  diluted  water  and  our  prediction  holds;  however,  at  higher 
bubble  concentrations  the  mixture  tends  to  rapidly  become  non-Newtonian, 
exhibiting  a  yield  stress,  and  a  decreasing  apparent  viscosity  with  increasing 
shear  rate.  Wallis  states  that  this  transition  occurs  at  void  fraction  above  about 
0.05  (Specific  Gravity  =  0.95).  We  believe  to  have  acquired  empirical  data  to 
support  the  theory  and  although  we  consider  our  findings  to  be  very  interesting, 
further  investigation  is  outside  the  scope  of  this  thesis  and  we  recommend  it  for 
future  work. 
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VI.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

In  our  experimental  investigation  we  show  that  the  introduction  of  gas 
bubbles  into  a  liquid  decreases  the  average  density,  and  thus  decreases  the 
buoyant  force  on  a  floating  body.  We  investigated  the  critical  average  density 
required  to  sink  a  buoyant  body  in  water  with  bubbles.  A  theory  of  the  critical 
density  for  sinking  was  developed,  and  predicted  that  the  average  fluid  density  is 
greater  than  the  ball  density  for  sinking.  The  experimental  data,  which  include  a 
quantitative  error  analysis,  agreed  well  with  the  theory  for  average  ball  densities 
at  low  flow  rates,  but  show  a  definite  trend  of  fluid  densities  that  are  smaller  than 
those  predicted  for  high  flow  rates.  This  may  be  due  to  the  abrupt  onset  of  non- 
Newtonian  behavior  of  the  fluid. 


B.  FUTURE  WORK 

This  research  would  benefit  from  substantial  further  study.  This  thesis 
investigated  the  bounded  case  where  bubbles  are  generated  uniformly  over  the 
cross-section  of  a  container.  Follow  on  research  could  investigate  the  case 
where  the  cross-sectional  area  of  the  liquid  is  much  greater  than  that  of  the 
bubbles,  which  would  occur  in  the  ocean,  for  example.  A  preliminary  qualitative 
experiment  showed  that  a  far  greater  flowrate  was  required  to  sink  the  body  in  a 
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more  open  environment  then  was  needed  during  our  experiment.  Also  the 
transitional  region  should  be  further  investigated  to  confirm  our  speculation  that  it 
is  indeed  a  change  from  Newtonian  to  non-Newtonian  fluid  characteristics. 
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APPENDIX  A.  DEMONSTRATIONS 


Our  investigations  of  the  introduction  of  bubbles  to  sink  a  body  floating  in  a 
liquid  naturally  began  with  ice  and  water.  We  were  partly  motivated  to  do  this  as 
a  dramatic  lecture  demonstration  that  would  impress  introductory  physics 
students.  We  employed  a  specially  designed  acrylic  cylinder  of  height  2.0  feet 
and  inner  diameter  3.0  inch  (Fig.  A.1 ).  The  cylinder  is  open  at  the  top  and  closed 
at  the  bottom,  and  has  an  air  inlet  below  an  acrylic  plate  that  is  4.0  inches  from 
the  bottom.  Drilled  through  this  plate  are  97  holes  of  diameter  1/32  inch,  and 
uniformly  spaced  1/4  inch  apart.  It  should  be  noted  that  the  air  chamber  beneath 
the  plate  can  be  pressurized  before  water  is  added  to  the  cylinder,  in  order  to 
prevent  the  water  from  leaking  through  the  holes.  An  alternative  is  to  close  the 
air  inlet  with  a  valve,  because  water  then  enters  the  chamber  at  a  sufficiently 
slow  rate  that  there  is  ample  time  to  perform  the  demonstration. 

Ice  cubes  roughly  1  inch  on  a  side  were  used.  To  be  visible  to  more  than 
a  small  audience,  the  ice  should  be  prepared  with  red  food  coloring  or  dye.  With 
the  apparatus,  we  are  able  to  sink  the  ice  with  sufficient  flow  from  a  compressed 
air  cylinder.  However,  after  the  ice  sinks,  it  continues  to  bob  up  and  down  due  to 
turbulence,  and  occasionally  even  returns  to  the  surface  for  a  short  time.  The 
demonstration  is  more  easily  and  convincingly  performed  with  a  body  whose 
average  density  is  very  near  that  of  water.  Such  a  body  sinks  at  a  lower  flow 
rate,  which  causes  less  turbulence.  In  addition,  the  air  can  be  simply  supplied 
from  the  demonstrator’s  lungs  to  sink  the  body.  A  small  closed  bottle  that  has 


37 


been  appropriately  weighted  works  well. 

Permanent  sinking  can  be  shown  by  using  a  body  similar  to  a  boat,  in 
which  the  material  of  the  body  has  an  average  density  greater  than  that  of  water. 
For  example,  a  metal  cap  from  a  bottle  can  be  employed.  However,  it  is  not  clear 
that  the  sinking  in  this  case  is  due  to  a  lowering  of  the  average  density  due  to  the 
bubbles,  because  the  substantial  foaming  that  occurs  on  the  surface  causes 
water  to  collect  inside  the  cap. 
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Figure  A.1 .  Original  apparatus  employed  to  sink  a  cube  of  ice  with  the  introduction  of  bubbles. 

Motivated  to  construct  an  apparatus  that  is  more  easily  observed  and 
controlled,  we  used  a  commercially-available  “density  ball,”  which  is  a  hollow, 
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stainless  steel,  10.2  cm  diameter  sphere  with  a  removable  plug  so  that  liquid  can 
be  added  or  removed  to  vary  the  ball’s  mass  and  thus  its  average  density  (Cenco 
Scientific  Co.  cat.  #76595N).  Because  the  ball  was  too  large  to  fit  in  the  above 
apparatus,  we  experimented  with  several  other  bubble-making  devices  in  larger 
containers.  These  “bubblers”  included  coiled  flexible  tubing  in  which  small  holes 
were  made  by  puncturing,  and  “frits”  (aggregates  of  small  clumps  of  metal).  The 
apparatus  that  provided  the  most  uniform  bubble  field  with  little  turbulence  and 
small  bubbles  was  a  set  of  four  aquarium  bubblers,  (Regent  PL-T714)  which  are 
hollow  porous  tubes  capped  at  one  end  and  with  an  air  inlet  at  the  other  (Fig. 

A.2).  We  cut  appropriate  lengths  to  span  the  bottom  of  a  4.0-liter  beaker.  The 
bubblers  are  held  on  the  bottom  of  the  beaker  with  suction  cups  that  clip  onto  the 
bubblers.  Air  is  supplied  by  a  compressed  air  cylinder  through  a  pressure 
regulator,  and  distributed  with  a  four-way  controllable  manifold  designed  for 
aquarium  use.  By  adjusting  the  average  density  of  the  ball  to  be  a  few  percent 
less  than  the  density  of  water,  we  are  able  to  sink  the  ball  with  a  flow  rate 
sufficiently  low  to  cause  little  turbulence.  When  the  density  of  the  ball  is  less  than 
about  95%  of  the  density  of  water,  flow  rates  high  enough  to  sink  the  ball  cause 
noticeable  turbulence.  This  apparatus  functions  well,  and  can  be  observed  even 
in  a  large  lecture  hall. 


39 


From  air 
source 


0 


APPENDIX  B.  ERROR  ANALYSIS  DETAILS 


I.  FLOWMETER  SPECIFICATIONS 


Table  B.1  provides  the  flowmeter  specifications  that  are  used  in  the  error 
analysis. 


Manufacturer 

Flow  Rate  Range 

Repeatability 

Resolution 

(liters/min) 

(liters/min) 

Gilmont 

0-16 

±  0.25%  Reading 

±0.25 

Gilmont 

0-44 

±  0.25%  Reading 

±0.25 

McMillian 

40  -  200 

±  1  %  Full  Scale 

±0.50 

Table  B.1.  Flowmel 

ter  Specifications. 

II.  COLUMN  DEFINITIONS 

The  following  is  an  explanation  of  the  various  columns  found  in  Table  B.2. 
The  resulting  data  is  used  in  the  error  analysis  to  determine  the  error  bars  as 
shown  in  Fig.  (5.4),  (5.5),  and  (5.6).  Upper-case  letters  designate  the  columns. 

A:  Specification  of  the  flowmeter  employed  for  the  current  trial. 

B:  Value  of  the  average  specific  gravity  to  which  the  density  ball  is 

adjusted. 

C:  Specific  gravity  below  which  the  ball  is  predicted  to  sink  (Eq.  (2.3)). 
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D:  Specific  gravity  at  which  the  ball  actually  sank  (Eq.  (4.8)). 

E:  Total  mass  of  ball  (including  liquid  inside). 

F:  Flow  rate  for  which  the  ball  sinks  below  the  surface  was  recorded. 

This  is  an  average  for  seven  separate  trails. 

G:  Flow  rate  for  which  the  ball  sinks  below  the  surface  was  recorded. 

This  is  the  maximum  of  seven  trials. 

H:  Flow  rate  for  which  the  ball  sinks  below  the  surface  was  recorded. 

This  is  the  minimum  of  seven  trails. 

J:  Flow  meter  repeatability  as  given  in  Table  B.1  at  the  average 

flowrate. 

K:  Flowmeter  resolution  as  given  in  Table  B.1  at  the  average  flowrate. 

L:  Flow  rate  for  which  the  ball  sinks  below  the  surface  was  recorded. 

This  is  the  half  of  the  difference  between  the  maximum  flow  rate 

G  —  H 

and  the  minimum  flow  rate - , 

2 

where  G,  and  H  are  the  results  form  columns  G,  and  H  respectively. 
M:  The  uncertainty  in  measuring  the  flow  rate  at  which  the  ball  actually 

sank.  VJ2  +K2  +  1}  , 

where  J,  K,  and  H  are  the  results  form  columns  J,  K,  and  L 
respectively. 
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N:  The  uncertainty  in  measuring  the  height  when  the  bubbles  are 

turned  on  and  steady  state  has  been  reached.  We  assess  this 
uncertainty  to  be  ±0.25  mm ,  independent  of  flowrate. 

O:  Flow  rates  ranging  from  0  to  1 39  l/min,  at  1 5  different  values,  were 

set  at  the  nitrogen  source,  and  the  corresponding  changes  in  height 
were  recorded.  A  3rd  order  polynomial  was  used  in  a  least  squares 
fit  through  the  points 

f  ( A/7)  =  -0. 1 71 507  +  0.41 9662(Q)  +  3.330407  x  1 0‘4  ( Q f 
+3.297877  x  1 0“5  (Q)3. 

This  was  evaluated  at  the  average  flow  rate  Q  =  Qavg 
P:  The  change  in  f(Ah)  with  respect  to  the  change  in  the  average 

flowrate  Qavg . 

=  9.89363 E  -  5 (Q)2  +  6.66081E  -  4(Q)  +  0.41 9662  . 
dQ 

This  was  evaluated  at  the  average  flow  rate  Q  =  Qavg 
Q:  The  uncertainty  in  the  change  in  height. 

6 (Ah)  =  yj(Shf  +  (df(Ah)/dQ)2(6Q)z . 

R:  The  overall  uncertainty  of  the  fluid  specific  gravity. 

Sp_  _  1  S  (Ah) 

pw  (1  +  Ah/h0f  h0 
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